SNUTP-94/13 (Revised) 
Heavy Quark Symmetry and the Skyrme Model 



Yongseok Oh 
Department of Physics, National Taiwan University, Taipei, Taiwan 10764, R-O.C. 



Byung-Yoon Park 
\^+-~ ' Department of Physics, Chungnam National University, Daejeon 305-764, Korea 

o\' 

Dong-Pil Min 

Center for Theoretical Physics and Department of Physics, Seoul National University, Seoul 

151-742, Korea 

>. 

^t ■ Abstract 

(N 

We present a consistent way of describing heavy baryons containing a heavy 
quark as bound states of an SU(2) soliton and heavy mesons. The resulting 
mass formula reveals the heavy quark symmetry explicitly. By extending 
the model to the orbitally excited states, we establish the generic structure 
of the heavy baryon spectrum. As anticipated from the heavy quark spin 
QJ \ symmetry, the c-factor denoting the hyperfine splitting constant vanishes and 

the baryons with the same angular momentum of light degrees of freedom 
form degenerate doublets. This approach is also applied to the pentaquark 
^ ■ exotic baryons, where the conventional c-factor plays no more a role of the 

hyperfine constant. After diagonalizing the Hamiltonian of order N^ 1 , we 
get the degenerate doublets, which implies the vanishing of genuine hyperfine 
splitting. 
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I. INTRODUCTION 

Hadrons containing a single heavy quark (Q) with its mass (tuq) much greater than a 
typical scale of strong interactions (Aqqq) can be viewed as a freely propagating point-like 
color source dressed by light degrees of freedom such as light quarks and gluons. Besides the 
chiral symmetry for the light quark system, such a system reveals an additional symmetry, 
so-called the heavy quark spin-flavor symmetry [l|-|3[ as the heavy quark mass goes to infinity. 
In this limit, the heavy quark spin decouples to the rest of the strongly interacting light quark 
system, since their coupling is a relativistic effect of order l/rug. Furthermore, since the 
heavy quark can hardly change its velocity due to the strong interactions via soft gluons, the 
dynamics of the system is independent of its mass and, therefore, its flavor. As a consequence 
of the heavy quark spin symmetry, the hadrons come in degenerate doublets M with total 
spin (regardless of the number of light quarks) 

3±=3i±h (1-1) 

(unless ji = 0) which are formed by combining the spin of the heavy quark with the total 
angular momentum of the light degrees of freedom ji- in other words, the total angular 
momentum of the light degrees of freedom, Ji = J — Sq, is conserved and the corresponding 
quantum number ji can classify the hadrons with a single heavy quark. 

In the Skyrme model a la Callan and Klebanov (CK) |J, heavy baryons can be described 
by bound states of a soliton of the 577(2) chiral Lagrangian and the heavy meson containing 
the heavy quark. This picture which was originally suggested for the strange hyperons has 
been shown to work successfully in describing the static properties of the heavy baryons 
with charm(c)- or bottom(b)-quark ||-^|. It is also extended to the study of the exotic 
pentaquark-baryons (P-baryons, in short) as the bound states of soliton and antiflavored 
heavy mesons 0. (Here, by "antiflavored" heavy mesons, we denote the Qq heavy mesons 
carrying opposite heavy flavor, i.e. , C = — 1 or B = +1. They are the antiparticles of 
normal heavy mesons consist of Qq, which we will call as heavy mesons.) However, one 
should be careful in applying the bound state approach to the heavy baryons, especially in 
case of P-baryons, due to following respects: 

Firstly, one should be more careful on the interpretation and the calculation of the so- 
called hyperfme constant: When quantized, the bound system of the soliton and heavy 
mesons with energy uob and grand spin quantum number k (the corresponding operator 
being defined as K — Ih + Jh', Ih'- isospin of heavy mesons, Jh- spin of heavy mesons) can be 
identified with the heavy baryon of isospin i and spin j. The masses of such heavy baryons 
in the conventional approach || turns out to be 

m (M . } = M sol + u B + ^{cj(j + !) + (!- c)i(i + 1) + c(c - l)k(k + 1)}, (1.2) 

where M so i and X are the soliton mass and its moment of inertia with respect to the collective 
isospin rotation, respectively. The spin of the heavy baryon j takes a priori one of the values 
\i — k\,- ■ ■ ,i + k and c, "hyperfme constant", is a constant defined through 

(k,k' 3 \e\k,k 3 ) = -c(k,k' 3 \K\k,k 3 ), (1.3) 



on the analogy of the Lande's g-factor in atomic physics. Here, G is the meson field operator 
induced by the collective rotation, which forms the first rank tensor in the space of the grand 
spin eigenstates \k,ks). The c yields the hyperfine splittings between the heavy baryon 
masses and it has been known to play the role of an order parameter for the heavy quark 
symmetry; it is required to vanish in the heavy meson mass limit so that the heavy baryon 
masses do not depend on the total spin. This requirement for the c, however, has been 
over- emphasized. In fact, the heavy quark symmetry does not necessarily require such an 
entire independence of mass on the total spin, as will be discussed in this work. 

To heavy baryons carrying a heavy quark one may extend straightforwardly the applica- 
bility of Eq.( |1.2| ). The heavy quark symmetry is manifested by vanishing c-factors of bound 
heavy meson states, provided a subtle point is corrected | lQfl . In the strangeness sector the 



vector mesons K* are eliminated out in favor of the pseudoscalar meson K, in analogy to 



the p mesons for the light quark system. [11] This approximation is valid only when the 



vector meson masses are sufficiently larger than those of the corresponding pseudoscalar 
mesons. As the heavy quark mass in both mesons increases, the heavy vector meson and 
the heavy pseudoscalar meson should be treated on the same footing. Defects in taking the 
conventional bound state approach of CK to heavy baryons have been pointed out and cor- 
rect heavy baryon mass spectra with the explicit heavy quark symmetry have been obtained 
by Jenkins et al. ]T2| and other groups []T^-[r6| in the infinitely heavy meson mass limit. 



In Refs. Hl2| , |13|| , a different but equivalent quantization scheme is adopted: the soliton is 
first quantized to nucleons and A's and then the heavy mesons are bound to them to form 
a heavy baryon, while in conventional bound state approach || the whole soliton-heavy 
meson bound system is quantized by using the collective coordinates. The mass predictions 
1~1| . [T5|] based on Eq.(|1.2|), however, give slightly different results from those of Ref. JT2J] , 



i.e., constant shifts of heavy baryon masses to the amount of 3/8X. This difference comes 
from the approximate treatment on the last term of Eq.( |1.2| ). We will show that the two 
predictions are equivalent by treating the term in proper way. 

We have met a completely different situation in the study of penta-quark heavy baryons 
which carry one heavy anti-quark. |17] The c-factors associated with the bound states of the 



antiflavored heavy mesons do not vanish. In some specific cases, one gets involved with a 
serious problem such as negative c-factors, with which the mass formula ( |1.2j ) would yield 
lower mass for the baryon with higher spin. However, this flaw can be cured by the following 
observation. There appear multiply degenerate heavy meson bound states, which prevent 
us from using the mass formula (|1.2|) in the present form. The c-factor defined through 
Eq.( |1.3| ) cannot play the role of the hyperfine constant. Therefore, the hyperfine splitting 
being of order 1/N C should be obtained, due to the appearance of the off-diagonal terms, by 
diagonalizing the Hamiltonian matrix with respect to the soliton-antiflavored heavy meson 
bound states degenerate up to order N®. We shall see that the mixing of the states at the 
N^ 1 order cannot be neglected in anti-flavor case. The resulting P-baryon masses do respect 
the heavy quark symmetry and the resulting hyperfine constants vanish. 

The main purpose of this paper is to clarify the things associated with the hyperfine 
constant in the conventional bound state approach and develop a consistent bound state 
approach to be applied not only to the normal heavy baryons but also to negative parity 
heavy baryons and exotic baryons carrying an antiflavor. It also supplies some of the details 



left out in our previous paper |L7j where we investigated the pentaquark exotic baryons. We 
will work in an extreme limit where both the soliton and heavy mesons are infinitely heavy 
and sit on the same point in space. This approximation enables us to get useful informations 
without getting involved with any complicated numerical calculations. 

This paper is organized as follows. In Sec. II, we briefly introduce our working Lagrangian 
density. The positive and negative parity eigenstates of the heavy mesons under the static 
potentials provided by the soliton configuration are found in Sec. III. Section IV is to discuss 
the collective coordinate quantization procedure for the bound system of a single heavy 
meson to the soliton. In Sec. V we describe the heavy baryons as the bound heavy mesons 
to soliton. We also derive a mass formula for the heavy baryons, which is more appropriate 
to appreciate the heavy quark symmetry than Eq. (|1.2|) . The realization of heavy quark 
symmetry in heavy baryon spectrum is also discussed. We study the pentaquark exotic 
baryons by considering the bound states of the "anti-flavored" heavy mesons in Sec. VI. We 
shall show that there exist degenerate doublets as given in Eq.flLlp in pentaquark states. 
A few concluding remarks are given in Sec. VII and explicit formulas are provided in 
Appendices A and B. 

II. HEAVY MESON EFFECTIVE LAGRANGIAN 

We start with describing briefly the effective Lagrangian for the heavy mesons interacting 
with Goldstone bosons, which respects both heavy quark symmetry and chiral SU(2)l x 
SU(2)$\ symmetry. (See Refs. 0HH for details.) 



Consider heavy mesons containing a heavy quark Q and a light antiquark q. Here, the 
light antiquark in a heavy meson is assumed to form a point-like object with the heavy 
quark, endowing it with appropriate color, flavor, angular momentum and parity. Let $ 
and $* be the field operators that annihilate j n =0~ and 1~ heavy mesons with C = +1 or 
B = — 1. These fields $ and $* form an 577(2) antidoublets: for example, when the heavy 
quark constituent is the c-quark, 



<5> = (D°,D + ) (& = {D°,D 



)-\T 



(2-1) 
$* = (_D*° 7 D* + ) ( $*t = (D*°, D*~ yj 



The traditional Lagrangian for the free fields is 

£ free = d^d»& - m|$$t - |$*^$*t, + m|*$*"$* t , (2.2) 

where $*^ = d^* u — ct$* is the field strength tensor of the heavy vector meson fields and 
m<j, and m$» are the masses of heavy pseudoscalar and vector mesons, respectively. 

In the limit of infinite heavy quark mass, the heavy quark symmetry implies that the 
dynamics of the heavy mesons depends trivially on their spin and mass. Such trivial depen- 
dence can be eliminated by introducing a redefined 4x4 matrix field H(x) as PJ2T 



l We will work with two light flavors. For the generalization to three flavors, see Ref. [21] 



H = ±±£ (^76 " K,r) , (2-3) 

instead of the traditional heavy meson fields, $ and $*. Here, we use the conventional Dirac 
7-matrices and j> denotes v fl r y fl . The fields $„ and <I>* , respectively, represent the heavy 
pseudoscalar field and heavy vector fields in the moving frame with a four velocity v^. They 
are related to the $ and $* as @ 



3» 



d> — r -iv-xm& 1 (f, 
<f>* __ p -iv-xm^* 1 (f,* 

Under the heavy quark spin rotation, H transforms 

H ^SH, 

with S G SU(2) V (the heavy quark spin symmetry group boosted by the velocity v). In the 
heavy meson rest frame, i.e. , t> M = (1, 0), S can be written explicitly as 

H -> e'^if, (p|c) 

where <Tj = §£yfc[7j, 7&], the Dirac spin matrices. To leading order in the heavy meson mass, 
the free Lagrangian density describing the heavy mesons propagating with a four-velocity 
v^ is nothing but that for the freely propagating heavy quark ||: 

£S* 8 = -ivpTtiPHH). (2.4) 

For later convenience, we have introduced H = 70^70, which transforms under the heavy 
quark spin rotation as H — > HS~ 1 . One may easily check that it comes as the leading order 
term in m$(= m$*) by substituting Eq. (|2.3| a) into £f rcc - 

On the other hand, the dynamics of the light quark system is governed by the SU(2)l x 
SU(2)r chiral symmetry, which is realized in a nonlinear way via a 2 x 2 unitary matrix 

, i ( tt V2ti + \ 
S = eXp( AUvr- -*°J)' (2 - 5) 

with the triplet of Goldstone bosons (tt + , tt° and tt~) and the pion decay constant /„.(=93 
MeV). Under the SU(2)l x SU(2)r transformation, £ transforms as 

E -> LZR\ (p|a) 

with global transformations L G SU(2)i and i? G SU(2)r. In terms of S, the interactions 
among the Goldstone bosons are described by the Lagrangian density 

C M =JfTr(d^d»Z) + ---, (2.6) 

where terms with more derivatives are abbreviated by the ellipsis. 

Construction of a chirally invariant Lagrangian for the couplings of the heavy mesons to 
the Goldstone bosons can be done by assigning H a proper transformation rule under the 



chiral transformation. There may be a considerable freedom. A standard one is to introduce 
a redefined matrix 



£ = VE, (2.7) 

which transforms under the SU(2) L x SU{2) R as 

£ _ L£U ] = U£R ] . (2.8) 

Here, U is a special unitary matrix depending on L, R and the Goldstone fields. From £ 
one can construct a vector field V^ and an axial vector field A^ 

which have simple chiral transformation properties: 

^ - UA^W. 



(2.9) 



(2.10) 

Let the light quark doublet q transforms as |24| 

q=(u,d) T -+Uq, (2.11) 

so that the heavy meson field H(x) transforms as 

H -»■ i/C/t. (2.12) 

This choice defines a simple transformation rule of H under the parity operation: 

£(r,*) - ^(-r,*), and F(r,*) -► 7 #(-f,£) 7 . (2.13) 

In terms of the vector field V^, a covariant derivative can be constructed as 

D,H(x) = H(d» +VJ), (2.14) 

which transforms under a chiral transformation as -D^if — > (D^H)U^ . 

Now, it is easy to write down a "heavy- quark- symmetric" and "chirally-invariant" La- 
grangian for the couplings of the heavy meson fields to the Goldstone bosons. To the leading 
order in the derivatives on the Goldstone boson fields, it reads 

Chqs = C M -tv,Ti(D^HH)+gTi(HYl 5 A^H), (2.15) 

with a universal coupling constant g for the $$*7r and $*$*7r interactions. The nonrela- 



tivistic quark model provides a naive estimation [19| for the value of g as 



9 = -I (2-16) 

On the other hand, the Lagrangian leads to the decay widths 

r(*;x - tt+$_i) = 2 r($;, -> n^ H ) = yj-f^l 3 , (2.17) 

J IT 

where the subscript ±| of $ and $* denotes the third component of their isospin. In case 



of Q=c, the c-quark, the experimental upper limit |25] of 131 keV on the D* width implies 
that |g| 2 ;$0.5 when combined with the D* + — > D + tt° and -D* + — > D°tt + branching ratios 
|26|| . In this work, however, we will take the pion decay constant f n and the heavy meson 



coupling constant g as adjustable parameters. 



III. SOLITON-HEAVY MESON BOUND STATE 

With a suitable stabilizing term,[] the nonlinear Lagrangian Cm supports a classical 



soliton solution 27 



£ (r) = exp(ir-fF(r)), (3.1) 

with the wavefunction F(r) satisfying the boundary conditions 

F(0) = vr and F(r) ™ 0. |J) 

The soliton solution carries a winding number identified with the baryon number 

(B. No.) = --j^ Jd 3 re l ' k Tr^E^-So^So) = 1, Qb) 

and a finite mass 

M so , = 4vr / o °°r 2 cir| (V 2 + 2^) + • ■ ■ , (0c) 

with F' = dF/dr and the ellipsis denoting the contributions from the soliton-stabilizing 
terms. 

Our main interest is the bound heavy-meson states (if any) due to the static potentials 
provided by the baryon-number-one soliton configuration ( |3.1| ) sitting at the origin. Note 
that we are working with the infinitely heavy soliton. Explicitly, the potentials are in the 
form of 

F = (y»F) = (0,4)fxf), 

A» = {A°,A) = (0,i(oi(r)r + o 2 (r)fT-f)), 



with 



, N sin 2 (F/2) , . sinF , . . . sinF 

v[r) = , ai(r) = and a 2 {r) = F 



In the rest frame f M = (1,0), the equation of motion for the eigenmodes H n (r)e~ t6nt of 
the iFfield can be read off as 

e n H n (r)=gH n (r)A-a } (3.3) 

where we have used that H(x)jj5 = —H(x)a. Here e n is the eigenenergy and n denotes 
a set of quantum numbers which classify the eigenmodes. The "hedgehog" configuration 
( ^.1[ ) correlates the isospin and the angular momentum, while the heavy-quark symmetry 
implies the heavy quark spin decoupling from the set. Thus, the equation of motion is 



2 The explicit form of the stabilizing terms is not essential for our discussions. Later we will adopt 
the Skyrme term as a stabilizing term for numerical results. For a review of the Skyrme model we 
refer to Ref. M. 



invariant under the parity operation, the heavy quark spin rotation, and the simultaneous 
rotations in the ensemble of spaces: isospin space, "light- quark spin" space, and ordinary 
space. Let L, St, Sq, and Ih be the orbital angular momentum, light quark spin, heavy 
quark spin, and isospin operators of the heavy mesons, respectively, and let Yi m (f), ^(±||, 
| ± |)q, and <fi±i be corresponding eigenstates, respectively (See Appendix A for their 
explicit representations.) The simultaneous rotations mentioned above are generated by the 
"light-quark grand spin" operators defined as 

K £ = L + S e + I h . (3.4) 

By the subscript £, we distinguish K% from the traditional grand spin operators (K = 
L + S + Ih with S = Se + Sq being the spin operators of the heavy mesons) used in the 
bound state approach || in the Skyrme model. Then, the eigenmodes of the heavy meson 
can be classified by the third component of heavy quark spin Sq, the grand spin {kg, k 3 ) and 
the parity it, so that n = \ki, k 3 , it, Sq}. 

The situation is very similar to obtaining the eigenmodes of the confined quarks in the 
chiral bag model ||29|| . We start with the construction of the eigenfunctions of the grand 
spin and the heavy quark spin by taking direct products of the four eigenstates, Ye me , (j>±k, 

£ (±||, and | ± \) Q , which yields four fC { k ] k , £=1,2,3,4: 



m 3 ,mt 



£&k3« =: H (£i,m e , l,m t \\i,m £ + m t )(Xi,mi + m t , ^,m s \ki, k3)Y ime (f)(p mt £ (m s \\s Q ) Q , 

(3.5) 



with the Clebsch-Gordan coefficients (£i, mi, £2, m<i\£,vn). Here, for a later convenience, 
we first combine the angular momentum and the isospin to form A(= L + Ih) and then 
combine the light quark spin with it. Although the heavy quark spin does not involve in the 
combination of the light-quark grand spin, we have included them in the definition of fC k % j k 8 

in order to shorten the expressions. The explicit forms of K, k l k are given in Appendix B. 

In terms of these K, k j k3S , the heavy meson wavefunction can be written as 

#n(r) = £ ^(r)41 v for vr = -(-1)^ states, 

i=1 ' 2 (3.6) 

Hn{?) = E h^{r)K,t} k for vr = +(-l) fe ^ states, 

t=3,4 Q 

with radial functions h\ . (r). Note that the electric modes (£=1,2 states) and the magnetic 
modes (i=3,4 states) are decoupled due to their different parities. 

Since the heavy mesons and the soliton are assumed infinitely heavy, their kinetic effects 
can be neglected and the heavy mesons are expected just to sit at the center of the soliton 
where the potentials have the lowest value. That is, in the heavy mass limit, all the radial 
functions h k ,(r) can be approximated as 

hl(r)= ai f(r), (3.7) 



with a constant «« and a function /(r) which is strongly peaked at the origin and normalized 
as J °° r 2 dr\f(r)\ 2 = 1. The problem is to find the eigenfunction of the equation 

eK k -, k3SQ = lgF'(0)JCkjk 3SQ {(r-r)[a-r\(T-r)}, (3.8) 

as a linear combination of ^tlk 3 s Q '-> l - e - > ^kjk 3 s a = Y,i a i^klk 3 s Q (^ = 1>2 or 3,4). In obtaining 
Eq.(|3.8|), we have used that F{r) = tt + F'(0)r + 0(r 3 ) near the origin so that ai(r) ~ 
-F'(0) + O(r 2 ), and a 2 (r) ~ 2F'(0) + O(r 2 ) and the identity (2a-rr-r - a-f) = (r-f ) (<t-t) (r-f) . 
The expansion coefficients a« are obtained by solving the secular equation 

Y J M ij a j = -ea h (i, j=l,2 or 3,4) (3.9) 

i 

where the matrix elements A4ij(i,j=l,2 or 3,4) are defined as 

M i3 = ^^0)/^Tr {4^ 3S Jr-f)[(a-f)](r-f)4f fc3SQ }, (3.10) 

with JC k J ek3S = 7°/C^ feas 7°. The minus sign in the right hand side of Eq.( |3.9|) comes from 



the fact that the basis states JC k k are normalized as (|B4| ) 



Explicit matrix elements are presented in Appendix B, according to which /Q *. s (f) are 

already the eigenstates of Eq. (|3.8|) with the degenerate eigenenergy ^gF'(0) and M. should 
be diagonalized for i, j = 3,4. The diagonalization of .M^, j=3,4) leads to two eigenstates: 

^k e k 3 s Q = V 2k^T^hk 3 s Q + V 2ft+l^k e k 3 s Q ' ( £ = +29 F '(0)) 
fck e k 3SQ = V 2k^+l^hk 3S Q - V 2k^l^hk 3 s Q - ( £ = ~2dF(0)) 

Thus, for each set of different quantum numbers {kii^f 0), k 3 , Sq}, f\ we have one state with 
the eigenenergy e = — |^F'(0) and three degenerate states with e = |gF'(0). Since g < 
and F'(0) < (in case of baryon- number-one soliton solution), we have one soliton-heavy 
meson bound state with a binding energy ^gF'(0). The positive energy states imply three 
(one for ki = 0) soliton- antiflavored heavy meson bound states with a binding energy ^gF'(0). 
|j~7| (See Sec. VI) In Fig. 1 given are energy levels of the bound states of heavy mesons and 
antiflavored ones. Shaded areas denote the continuum states and triply degenerate states 
are represented as thick lines. 



In case of a typical soliton solution stabilized by the Skyrme term [p7], we have -^'(0) 



-2ef v with e being the Skyrme parameter. When the parameters are fixed as / 7r =64.5 MeV 



3 In case of ki = 0, we have two eigenstates; 



4oi Q (r): s = +yF'(0), 
4ol Q (r): s = -lgF'(0). 



and e=5.45 for the soliton to fit well the nucleoli and Delta masses ||30| , F'(0) amounts 



to ~ —0.70 GeV which implies that the binding energies of the soliton-heavy meson and 
soliton-antiflavored heavy meson bound states are ~gF'(0) ~ 0.79 GeV and |gF'(0) ~ 0.26 
GeV, respectively. Compared with those of Refs. |||9j , one can see that the binding energies 
are reduced by a factor half and more. It should be emphasized further that in Refs. |||9| 
the binding energy increases as the heavy meson mass increases and also that our results 
are obtained with infinite heavy meson mass. 

Each state can be combined with the heavy quark spin to produce doubly-degenerate 
grand spin eigenstates with k± = ki ± 1/2 (provided kg ^ 0, for which we have only a grand 
spin k = 1/2 state). These degeneracies are the direct consequence of the heavy quark 
symmetry and a special attention should be paid to the quantization procedure, although 
the degeneracy in kg (and in principal quantum numbers which are related with the radial 
excitations) may be an artifact from the approximation ( |3.7| ) on the radial function hjl ' (r). 
In general, when the heavy meson's kinetic term is taken into account, the radial function 
feels the centrifugal potential £ e s(£ e s + l)/r 2 near the origin so that it behaves as h^' ~ r e cs. 
Here, £ e ff is the "effective" angular momentum ||, which is related to £ as 



Kcff 



+ 1, if A = £ + i 
- 1, if A = i 



I 1 (3.12) 

2' 



Due to the vector potential V (~ i(f x r)/r, near the origin) the singular structure of 
D 2 = (V - Vf is altered from £{£ + l)/r 2 (that of V 2 ) to £ eS (£ eS + l)/r 2 . Thus, only 
those states with £ e g = can have strongly peaked radial function and the degeneracies will 
be broken in such a way that the states with higher £ e s have the higher energy. For the 
positive parity states, £ e g = can be achieved only when £ = 1. In Table I listed are a few 
heavy meson eigenstates which involve £=0,1 angular basis. From now on, we will restrict 
our considerations to these states. Note that the £ e ^ of the wavefunctions /C^ s , /QnL 

(included in K} 2 rns )> an d K-2m,s (included in /Q ms ) is 2 and that the /C^ s (included in 

ICims ) is 1. In case of the finite heavy meson mass, these states will have higher eigenenergies 
and will become even unbound. 

In terms of the eigenmodes, we can expand the heavy meson field operator as 

H(x) = J2H n (^e~ l£ - t a nj (3.13) 

n 

with the heavy meson annihilation operators a n . Note that we don't need to include the term 
for the antiparticles. The Fock states on which the quark field operators act are obtained as 

|m,n 2 ,---) =4 1 4 2 ---|vac), (3.14) 

where |vac) is the vacuum of the heavy meson fields. Hereafter, we will denote the Fock 
states of a single heavy meson occupying the corresponding state as in Table I. (To simplify 
the notations, unless necessary, we will not specify such trivial quantum numbers as the 
third component of the grand spin, the heavy quark spin and the parity; k%, sq and it.) 



10 



IV. COLLECTIVE COORDINATE QUANTIZATION 

What we have obtained so far is the soliton-heavy meson (or antiflavored heavy meson) 
bound state which carries a baryon number and a heavy flavor (or anti-heavy flavor). To 
endow the states with correct quantum numbers such as spin and isospin, we have to go to 
next order in 1/N C , while remaining in the same order in ttiq, namely C^m-giV" 1 ). This can 
be done by quantizing the zero modes associated with the invariance under simultaneous 
SU(2) rotation of the soliton configuration together with the heavy meson fields: 

So-tC&Cft, and H^HC\ (4.1) 

with an arbitrary constant £77(2) matrix C and £q = Uq. The rotation becomes dynamical 
by giving time dependence to the £77(2) collective variables as 

Z(?,t) = C(t)W> f )&(t), and H(r,t) = H bi (r,t)<y(t), (4.2) 

and then the quantization is done by elevating the collective variables to the corresponding 
quantum mechanical operators. In Eq.( |4.2j ), H^ refers to the heavy meson field in the 
isospin-co-moving frame, while H(r, t) refers to that in the laboratory frame. Substitution 
of Eq. (|4.2j ) into Eq. (|2.15|) leads us to the Lagrangian (in the reference frame where the heavy 



meson is at rest in space but rotating in isospin space) 

Z/ ot = -M sol + fd 3 r {-iTr(d H h{ H h{ ) + gTr(H M A-a H M )\ 

f r , ( 4 - 3 ) 

+±Zcu 2 -ljd 3 r{ Tr(77 bf |(£tf . ££ + £f . J?) H M )} , 

where we have kept terms up to 0{w^qN~ 1 ). The "angular velocity", uj, of the collective 
rotation is defined by 

CftdoC = \it-Q } (|Ja) 

and X is the moment of inertia of the soliton configuration with respect to the rotation 

X = f fl j r 2 dr(sm 2 F + ■ ■ ■) §Jb) 

where the contributions from the soliton-stabilizing-Lagrangian are abbreviated simply by 
ellipsis. 

Given the Lagrangian ([Of) that describes dynamics up to order 0(mQN~ l ), one has 
equation of motion 

id H u = H M [gA-a - \^r ■ ^ + & • $&)] (4-4) 

consistent to that order. The last "Coriolis" term in the equation of motion couples the fast 
and slow degrees of freedom. Although the heavy mesons are infinitely heavy, their angular 
momentum and isospin are associated with the light constituents. Thus, we may take those 
light degrees of freedom of the heavy meson fields as the fast ones and the collective rotations 
as the slow ones. Note that the scale of the eigenenergies \e n \ of the heavy mesons is much 
greater than that of the rotational velocity; \e n \ 3> |a>|. 
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Generally accepted procedure of handling these different scale is as follows. Jl5| We 
first solve the equation of motion for fast degrees of freedom with slow degrees of freedom 
"frozen" . In this way, we get "snap-shot" pictures of the fast motion. Next we solve the 
equation of motion for slow degrees of freedom taking into account the "relic" of the fast 
motion that has been "integrated out", in a manner completely analogous to the incor- 



poration of Berry phases [pi] . It is also analogous to the "strong-coupling limit" of the 
particle-rotor model |£| in nuclear physics, where the coupling between the rotating "core" 
and the particle is much stronger than the perturbation of the single-particle motion by 
Coriolis interaction. Here, the roles of the particle and the rotor are played by the bound 
heavy-mesons and the rotating soliton configuration. Thus, we may take the assumption 
that the bound heavy mesons rotate together with the soliton core in the unchanged eigen- 
modes. It enables us to expand the H^(x) in terms of the classical eigenmodes obtained in 
Sec. Ill as Eg. (gig): 

H hi {x) = Y,H n {r)e-^ t a n . (4.5) 

n 

Taking Legendre transformation of the Lagrangian, we obtain the Hamiltonian as 
H= / dr \ 7777 7^bf, Q /3 + ~ x 0J a - L mt 



\S(H U ,ap) M ' Q/ 7 5uJ a (4.6) 

= M sol - gjd 3 r Tr(H hi A-aH M ) + ^(R - 6(oo)) 2 , 

— * 

where the rotor spin R is the canonical momenta conjugate to the collective variables C(t): 

5L TOt 
R a = -— = Tuo a + e o (oo), Ha) 

00J a 

with 6(oo) defined as 

0(oo) = -lJd 3 rTr[H M l^^ + ^)H h{ }, (0b) 

whose expectation value with respect to the state \n) is the Berry phase associated with 
the collective rotation [l4]|n| . Note that it is nothing but the isospin of the heavy mesons 
modulo the sign. (See Eq.( |A13| b) below) 

With the collective variable introduced in Eq.([4.2|), the isospin of the fields U(x) and 
H (x) is entirely shifted to C(t). To see this, consider the isospin rotation 

£ -► AY,A\ H -► HA\ (4.7) 

with A G SU(2)y, under which the collective variables and fields in body-fixed frame trans- 
form as 

C(t) -► AC(t), H hi (x) -. H hi (x). (4.8) 

Note that the if -field becomes isospin blind in the (isospin) co-moving frame. The conven- 
tional Noether construction leads the isospin of the system as 
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I* = \ Tr(r a Cr b tf){Iuj b + 6 (oo)) = D ab (C)R b , (4.9) 

with D ab (C) being the SU(2) adjoint representation associated with the collective variables 
C(t). ' 

Under a spatial rotation (together with the spin rotation in case of the heavy meson), 
with the help of the i^-symmetry, the fields transform as 

£(f,i) -> £(rV) = C(t)B^ (r)BC\t), 

H{f, t) -> e^H^r*, t)e"^^ = (e^^i^f*, t)e~ l ^ + ^)B^C(t), 

with r* = exp(«<3 • L)r and £> = exp(i|<3 • r) G SU{2). This means that the spatial rotation 
is equivalent to the transformation of the collective variables and if -fields in the body fixed 
frame as 

C(t) -* C(t)B\ , s 

i , (4.11) 

H hi (r,t) -> e i 2^ii bf (^,t)e- i 2°-(' ? + T? ). 

Therefore, we get the fact that the spin of the Hm(x) is the grand spin; that is, the isospin of 
the if -field is transmuted into the part of the spin in the isospin co-moving frame. Remember 
that the H^(x) becomes isospin blind in that frame. Applying the Noether theorem to the 
Lagrangian ( |4.3| ), we obtain the spin of the system explicitly as 

J = R + K hU (4.12) 

with the grand spin of the heavy meson fields in the isospin co-moving frame 

K hi = _ jtf r Tr{ (^ bf + [i^ Hh{] + # bf (_If))# bf }. (4.13) 

Finally, the heavy-quark spin symmetry of the Lagrangian under the transformation 

H{x) -► e l ^H{x) = {e l ^H M (x))C(t), (4.14) 

has nothing to do with the collective rotations. The heavy-quark spin operator remains 
unchanged in the isospin co-moving frame : 

S Q = - Jd 3 r Ti{\aHH) = - Jd 3 r Ti{\aH u H u ). (4.15) 

Because of this heavy-quark spin decoupling, it is convenient to proceed with the spin 
operators Je for the light degrees of freedom in the soliton-heavy meson bound system 
defined as 

J i = J-S Q = R + K e . (4.16) 

Upon canonical quantization, the collective variables become the quantum mechanical 
operators; the isospin (i), the spin (Jg) and the spin of the rotor (R) discussed so far become 
the corresponding operators i a , Jg ya and R a , respectively. We distinguish those operators 
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associated with the collective coordinate quantization by using a tilde on them. Let the 
eigenstates of the rotor-spin operator R a be denoted by |i;m 1 ,m 2 } (mi,m 2 = —i,—i + 
l,---,z): 

R 2 \i;mi,m 2 } = i(i + l)\i;m 1 ,m 2 }, 

R z \i;m 1 ,m 2 } = m 2 \r,m 1 ,m 2 }, (4-17) 

I z \i;m 1 ,m 2 } = mi\i;m 1 ,m 2 }. 

These states are represented explicitly by the Wigner D-functions 

V2i + lD® ;m2 (C). (4.18) 

At this point, it should be mentioned that we are following exactly the same quantization 
procedure of CK ||. We quantize the whole soliton-heavy meson bound system to obtain 



the heavy baryon states. In Refs. [|T2] , [I3|] , only the soliton is quantized to nucleons and A's 
by using the collective coordinate quantization. Then the heavy mesons with good isospin 
and spin are bound to form a heavy baryon. The corresponding isospin and spin operator 
of the heavy baryon system are different from Eqs. (|4.9| ) and (|4.12|) , respectively. However, 
both approaches lead to the same final results for the heavy baryon spectrum. 

V. HEAVY BARYONS 

The eigenstates \i,i3',je,jt } 3', Sq)) of the operators I a and J^ a with their corresponding 
quantum numbers i,i 3 (isospin) and jt,ji t 3 (spin of the light degrees of freedom) are given 
by the linear combinations of the direct product of the the rotor-spin eigenstates \i\ m 1 , m 2 } 
and the single-particle Fock state \n): 

\hh;jt,je,3;s Q )) a = J2(hJi,3-m } k^,m\j e} j i3 )\i;i 3 ,j i3 -m}\k e ,m,s Q } a . (5.1) 

m 

One may combine further the heavy quark spin and the spin of the light degrees of freedom 
to construct the states with a good total spin, which is, however, not necessary in our 
discussion. Remember that, in the infinite heavy quark mass limit, (je,ji,3) themselves are 
good quantum numbers to label the heavy hadrons. For a given set of (i,jt), there can be 
more than one state depending on which Fock state \n) is involved in the combination ( p.l| ). 
We will distinguish them by using a sequential number, a(=l,2- • •), in \i,ji))\. Here again, 
to shorten the expressions, we will not specify the quantum numbers i$, j$ and Sq unless 
necessary. 

In Table II, we list a few \i,jf)) states for low-lying heavy baryons. Here, we consider the 
\i,Je)) states made of the integer rotor-spin-states so that they describe the heavy baryons 
of a half-integer spin (J = ji ± |). In the Table, we do not present the states such as |0, 1 + )), 
|1,0 + )) and |0,0 - )), which cannot form the bound states. The incorporation of the 1/N C 
order corrections due to collective rotations cannot turn them into bound systems. We also 
exclude two other possible |1, 1~)) states from the list, since they come from the heavy-meson 
states |2) 12 of £ = 2 and thus they do not mix up with the state |1, 1~)) shown in Table II 
by the collective rotation. 
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The physical heavy baryons of our concern appear as the eigenstates of the Hamiltonian 
H . We will not try to find the exact eigenstates of the Hamiltonian. Remind that we have 
kept the terms only up to 0{m [ QN~ 1 ) in the collective coordinate quantization procedure. 
We take the last term in the Hamiltonian 

H rot =-L(R-Q(oo)) 2 , (5.2) 

as a perturbation of order 0(7tiqN~ 1 ) and we will consistently search for the approximate 
eigenstates to that order. 

Except for the case of |1, 1 + )), we have only one bound state for a given (i,jf). To the 
first order, the eigenstate of H is approximated by the unperturbed one \i,j£))i and the 
mass correction to the corresponding baryon is obtained by taking the expectation value of 
the H rot with respect to it: 



1 



Mr* = ((iJJ\(R- e(oo)) 2 |*,j7». (5.3) 



Now, our problem is reduced to evaluating the expectation values of the operators — 2R ■ 
0(oo) and 9 2 (oo). Using Eq.([B9l), we easily obtain the expectation value of the operator 
-2R-Q(oo) as 

a ((i,jJ\ - 1R ■ 0(oo)|i,tf» a = c aa {j e ( Je + 1) - i(i + 1) - h{h + 1)}, (5.4) 

with c aa being the c- value associated with the single-particle Fock state \kg) a participating 
in the construction of the state \i,jg))- As for the expectation value of the operator 2 (oo), 
in many CK-based models it has been approximated as 

((», jT|0 2 (oo)|z, jj)) = a {h\Q\^)\k t ) a « | a <k|3(oo)|A*> B | 2 = c 2 a Mh + 1). (5.5) 

In the heavy meson mass limit, it can be exactly evaluated as 

a (ke,k 3 \Q 2 (^)\ke,h) a = ]T | a (^,fc 3 |6(oo)|^,^) 6 | 2 , (5.6) 

{k' v k' 3 ,b} 

where the summation runs over the complete set of intermediate Fock states \k' e ,k' 3 ) b . In- 
cluding all the Fock states that have non-vanishing expectation value (n|0(oo)|m) with the 
help of the approximation ( J3.7] ) on the radial functions, we obtain 

a (^|0 2 (oo)|^) a = f. (5.7) 

One may obtain the same result by using the fact that 6(oo) defined as in Eq. (|4.6[ b) is 
nothing but the isospin operator Ih (modulo opposite sign) of the heavy meson field. 

Thus, for the heavy baryons with quantum numbers {i,jj) which allow only one bound 
\h3e)) state, we obtain the mass formula as 

3 1 
m {i , JV = M sol + uj b + — + —{cj £ (j e + 1) + (1 - c)i(i + 1) - ck e (k e + 1)}, (5.8) 
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where c is an abbreviation for c aa and ujb = m® — |g-F'(0). In order to compare it with 
Eq.( |1.2|) , we have included the weight averaged heavy-meson mass m$(= i(3m,j,* + m$)). 



Eq.( [5.8| ) is in a quite different form from Eq.( |1.2| ) and satisfies the heavy-quark symmetry 
regardless of the c-value. However, all the c-values associated with the heavy-meson bound 
states \ki)- vanish identically; viz., 

c ~ 2^TT C33 ~ 2 2k t + i ° 34 + 2^TT C44 " °' (5 ' 9) 



as given in Appendix B. Consequently, both mass formulas ( |1.2|) and ( |5.8| ) yield the same 
heavy baryon masses apart from a constant shift coming from the way of evaluating the 
expectation value of 2 (oo) as in ( |5.5| ) and ( |5.7| ). It should be emphasized that such a 
coincidence owes entirely to the vanishing c-values and the presence of the unique bound 
state. When one has no n- vanishing c-value as we shall see in Sec. VI, Eq. fll.2p cannot be 
applied anymore. 

If we have multiple degenerate bound states, \i,j][)) a (a = 1, 2, ■ • •), the situation becomes 
a little bit complicated. The mass corrections and the corresponding eigenstate are obtained 
by diagonalizing the energy matrix £ whose matrix element is defined by 

S ab = AhJl\H r ot\iJl)) b - (0,6 = 1,2, ■ ■ ■) (5.10) 

For the (i = l,jf = 1 + ) heavy baryons, we have doubly degenerate bound states; |l,l))i 
and |1,1)) 2 . Since they are, respectively, made of \kf = 0) 3 and \ki = 2)_, the first rank 
tensor 0(oo) cannot lead to nonvanishing energy matrix between the two states. Thus, each 
state can be separately the eigenstate with their degenerate massQ given by the mass formula 



Because of vanishing c-value, we can write the mass formula for the heavy baryon simply 
as 

1 3^ 



m (ijJ) = M ™i + "B + ^(i(i + 1) + f )• (5.11) 



In Fig. 2, the resulting heavy baryon spectrum is presented schematically. It is interesting 
to observe the degeneracy in mass of the heavy baryons with positive and negative parities. 
Since we are working upon an assumption on the radial functions with ignoring any effects 
of the kinetic term, we are not at a position to conclude whether such a parity doubling has 
any physical importance or it is just an artifact of the approximation. 
Explicitly, we have the masses of (0, + ) and (1, 1 + ) heavy baryons as 

771(0,0+) = M soi + uj b + 3/8J, 

(5.12) 

777.(1,1+) = M soi + u B + 11/8 J, 



4 The exact degeneracy is however an artifact of the approximation of using the same radial 
function for all the state. When the heavy mesons are allowed to move, 1 1 , 1 }} 2 will have higher 
energy. 
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which correspond to the mass of Ag and degenerate mass of Eq and Eg, respectively and are 
consistent with the results of Ref. |12|. Eq.( 5.12| ) yields an interesting model-independent 



relation for the mass difference of Eq(Eq) and Aq: 

m s Q ,s^ - m A Q = |( m A - m N ) « 0.20 GeV, (5.13) 

where we have used the fact that the SU(2) collective quantization of the bare soliton leads 
to the the nucleon and delta masses as 

3 15 

m N = M sol + — , and m A = M sol + —. (5.14) 

ox ox 

In the mass formula for the heavy baryons, we have three parameters, M so i, 1/2" and 
gF'(0). For a naive prediction on the heavy baryon masses, we fit them to produce ex- 
perimental values of ^.^(=939 MeV), m A (=1232 MeV) and m Ac (=2285 MeV), which leads 
to 

M so , = 866MeV, 1/X = 195 MeV and gF'(0) = 419 MeV. (5.15) 

Combined with the slope of the soliton wavefunction F'(0) ~ —690 MeV (in case of the 
Skyrme term-stabilized soliton solution), Eq. fl5.15|) implies g- value as g m —0.61 which 
is comparable to that of the non-relativistic quark model (—0.75) and the experimental 
estimation via the D*-decay (|g| 2 ;$0.5). This set of parameters yields a prediction on the 
the A b mass and the average mass of the E c -E* multiplets, ms c (= |(2m s . + rn^J) as 

m Ab = M sol + m B - fsF'(O) + 3/8X = 5623 MeV, 

(5.16) 

m Sc = M sol + m D - \gF'(ti) + 11/8 J = 2483MeV, 

which are comparable with the experimental value of the A^ mass (5641 MeV) and E c mass 
(2453 MeV) @. Recent experimental data for the E* mass of 2530 MeV @ (although it 
needs the verifications of other groups) gives 2504 MeV for the experimental value of m^, 
which is not far from our estimation. In our approach, the E c and E* are degenerate in 
mass. To get the splitting between them, one needs to include l/mq corrections Jl0| . 



With 1/X ~ 200 MeV, we can estimate the discrepancy in the masses of the heavy 
baryons given by the mass formulas, (|1.2j) and (|5.8|) , as 3/8X ~ 60 MeV. It amounts 10% 



of the binding energy —^gF'(0) ~ 630 MeV and 30% of the rotational energy ~ 1/X ~ 200 
MeV. One should not say that the discrepancy is at most 3% for the case of the charmed 
baryons by comparing it with the whole heavy baryon masses. Although we are working 
with the heavy baryons, our scheme is valid only in the low energy region below Aqqt> 

So far, we have considered the first order perturbations in the masses of the bound \i, j^)) 
states. To estimate naively the effects of the other states, we take into account the positive 
energy states into our procedure. It leads to a mass correction to the heavy baryon described 
by \iJt)) a as 

Am (y ,) = J2 la((Z,j7l ^° t|Z,jI))fc|2 + • • • , (5.17) 
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where the summation runs over the unbound states \i,Ji))b- They are, thus, at most of 
second order in 1/N C , which is out of our concern. Here, to check a consistency, we evaluate 
the leading order correction to (i = l,jj = 0~) heavy baryon state. We have 



^((i, o-|^ rot |i,o-)) 
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A ™d,o-) = - " ' ;;' — = -7tt^tt^7 ~ - 2 - 6 MeV > ( 5 - 18 ) 



2gF'(0) 2gF'(0) 18J 2 



which is negligibly small compared to the first order corrections, 11/8X ~ 270 MeV. The 
coupling of 1/N C order due to collective rotations cannot compete with the energy difference 
2gF'(0) of order iV c °. 

VI. PENTAQUARK EXOTIC BARYONS 

In the limit of infinite heavy quark mass, the quark model predicts stable pentaquark(P) 
exotic baryons whose quark contents are Qq 4 |35|j36|] . With the exact SU(3)p symmetry 
assumed for the light quarks, it was shown that a strange anti-charmed baryon P Ss (csq^, 
go = u, d) is stable against the decays into AD or ND S with binding energy about 150 MeV. 
The binding energy becomes down to ~ 85 MeV if included a realistic SU(3)f symmetry 



breaking |37j an d it becomes even unbound when the motions of the heavy constituent are 
taken into account P^pE0|. In the Skyrme model, as the heavy meson masses increases, 
there appears bound state (s) for the antiflavored heavy mesons to the soliton, which reveals 
a possibility for the stable nonstrange P-baryon(s) ||. It is interesting to note that in quark 
model such a nonstrange anticharmed baryon (cq^) cannot have sufficient symmetry to yield 
a hyperfine binding. 

Our approach can be easily switched to the one for the soliton-antiflavored heavy meson 
bound system by considering the negative energy solutions with the four velocity t> M in 
the equations replaced by — v^. (One may develop an effective Lagrangian proper for the 
antiflavored heavy mesons. See Ref. |T7]].) Now, for each k e (^ 0), we have three degenerate 
bound states of the antiflavored heavy mesons, \ke)i t 2,+, with the binding energy |gP'(0). 
Compared with that for the heavy mesons, the binding energy is reduced by a factor 3. 
With gF'(0) = 419 MeV as given by Eq.( |5.15| ), it amounts to ~ 210 MeV and is comparable 
to 1/X. 

In Table III listed are the bound \i,jj)) states proper for the discussions of the P-baryons. 
As for the states with quantum numbers (i = 0,jf = 0~), (1, 0~) and (0, 1~), we have only 
one bound state. To first order in 1/N C , the masses of such P-baryons are given by the same 
mass formula as Eq.( 



m( id7) = M sol +u' B + —{cjt(ji + 1) + (1 - c)i{i + 1) - ck e {k e + 1) + f }, (6.1) 

with oj' b = rrlq, — |pP'(0). The c-factors associated with |0)i and |1) + states are obtained 
as and —1/4, respectively. [f|l]] The negative c- value is remarkable. However, these states 
are associated with the zero rotor-spin state and such nonvanishing c-values do not play any 
important role in their masses which are simply obtained as 
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m 



(o,o-) — '"(o,i-) 



m (o,i-) = M soi + uj' b + —. (6.2) 



8J 



For the case of (1, ), the mass formula ( |6.1|) gives 

m ( i, -) = M sol +Lu' B + — = M N + m*- -gF'(0) + — . (6.3) 

Since ^gF'(0) ~ 210 MeV and 1/J ~ 195 MeV, the rotational energy blows up the state 
above the decay threshold of the nucleon-heavy meson bound system, M t h{= M N + m$). 

In cases of the (0, 1 + ) and (1, + ), we have two degenerate states |0, l + ))i,2 and |1, + ))i i2 , 
respectively. Since the states come from the single-particle Fock states of the same k^, the 
energy matrices can be expressed in a form of Eq. flb.lp with the constant c replaced by the 
2x2 matrix associated with kf — 1: 



i A 

Explicitly, we have 



°=li-2- < 6 - 4 ' 



Q 1 1 O 

£(o,i+) = M soi + u' B + — , and £ {lfi+) = M sol + J B + — - —. (6.5) 

Here again, the energy matrix ^(o,i+) is independent of c because the states are made of the 
zero rotor-spin state. Anyway, both are diagonal so that |0, l + ))i,2 and |l,0 + ))i ) 2 are the 
eigenstates of the Hamiltonian with the masses 

3 1 

m (o,i+)i, 2 = M soi + uj' B + — = M N + m$- -gF'(0), (6.6) 

and 

m (10+)l = M sol + u y B + — = M N +m 9 - -gF'(0) + — , 

19 12 (• ' 

m (1)0+)a = M sol + UB + — = M N + mz- -gF'(0) + -. 

The state |1,0)) 2 lies above the decay threshold. 

For the (i = l,jj = 1 + ) states, we have three possible combinations; |1, l + ))i,2,3- As far 
as the first two states |1, l + ))i,2 are concerned, the energy matrix reads simply 

%i+) = M Ml + a/ B + ^-|, (6.8) 



with the same 2x2 matrix c as in Eq.( |6.4| ). It leads us to the mass eigenvalues 

9 

m (i,i+)i = M *°i + ^ + 8J' 



(6.9) 



15 

m (l,l + )2 = M S0l+UJ B +—. 

With including the third state |1, l + ))3, the full 3x3 energy matrix is obtained as 

11 1 f" 10 ^ 
£(i,i+) = M sol + u' B + — + —\ 2 \. (6.10) 

81 4J V30 1 ) 
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The coupling between |1, l + ))i and |1, l + ))3 modifies m/ x 1+ n as 

7 

m (l,l+) = Ms ° l + U B + gj> ( 6 - n ) 

with |1, 1 + ))- = 2 |1> l + ))i — 1 1^-' 1 + )) 3 an< ^ adds doubly degenerate unbound states of mass 
M so i + u' B + 15/8X. Note that in obtaining those off-diagonal matrix elements one cannot 
use Eq.([B9|) and they cannot be written in a form of Eq.( |6.1| ) anymore. 

As for the |1,1~)) states, we have four possible combinations. They are made of the 
heavy-meson bound state with different kg(=0,l,2). By using Eq.([B7|), we obtain the energy 
matrix as 



% 1 -) = M s . + ^ + H + J_ 



/ 4v/3 \ 
A^/3 3 oVTE 

0-6 
\ v^O 9 J 



(6.12) 



The couplings between the states are somewhat strong. It yields four mass eigenenergies 

m (i,i-) = M sol + u B + — , 



i 
m 



M sol +u' B + ^ (6.13) 



(i,i-) ~~ soL ^ B "T" g J' 



2 
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m (Ci-) = Mso1 + u 'b + gj> 

where the states with m^ x _-j are doubly degenerate and expected to be bound. 

Our remarkably simple results on the P-baryon masses are summarized in Table IV 
and Fig. 3. To provide a rough scale, in Fig. 3, the heavy baryon spectrum is presented 
in the left hand side. Here again, the degeneracy of the P-baryon states in the parity is 
apparent. In Table IV, we also give a rough prediction obtained by using the values for 
the parameters given by Eq. (|5.15|) . However, all the states listed in Table IV do not seem 
to survive under the finite heavy meson mass corrections. Recently, we have reported that 
such finite mass corrections reduce the binding energy by an amount from 25% (in case of 
bottomed baryons) to 35% (in case of the charmed baryons) of their infinite mass limit, 
fgP'(O) |1(|. Note that 35% of fflP^O) is comparable to the binding energy \gF'{Q) for the 
soliton-antiflavored heavy mesons. The finite mass corrections will be more crucial for the 
heavy meson eigenstates with £ e s ^ 0. All the degeneracies in the heavy meson eigenstates 
will be broken. Then, the couplings between the states as represented by the energy matrices 
£(io+) j £(i,i-) an d £(i,i+) wm d° l ess important roles. However, it does not mean that one 
can apply the mass formula Eq. (|1.2j ) for the P-baryons in the present form. 

In order to show it, here, we compare our results with what would have been obtained by 
straightforwardly extending the bound state approach of CK as in Ref. 0. As a nontrivial 
example, we consider the P-baryons with i — 1 and j n = | , | (jj = + , 1 + ), which are 
obtained by combining the rotor-spin state with i = 1 to the heavy meson bound states with 
kj = 1 + . To simplify the process, we will include only the |l)i state into the consideration, 
in which case the P-baryon masses are simply given by the diagonal elements of the energy 
matrices £(i,o+) or £(i,i+) associated with |l,0 + ))i or |1, l + ))i : 
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m (i 1+) = M sol + u' B + (11 - 16ci)/8T, 
{ ' 2 ' (6.14) 

m (i.i+) = m (i,a+) = M soi + u'b + (11 _ 8ci)/8J, 

with ci = 4, the c- value associated with the state |l)i. In the traditional bound state 
approach with sufficiently heavy meson masses, 1 1} 2 will appear as nearly doubly degenerate 
states, say, \k = |) and \k = |), which become to completely degenerate in the infinitely 
heavy mass limit. When combined with i=l rotor spin state, the former yields P-baryons 
°f j v — 9 > I with masses 



m n 1+ =M sol +u' B + (11- 8c x )/8T. 

I 1 . 2 ) 2 



m 



(6.15) 

n 3+ =M soi + ^ + (ll + 4 Cl )/8X, 



and the latter yields 



m n 1+ . = M so , + u' B + (11 - 20c 3 )/8J, 
(.1,5 J 2 

m 3+ = M so , + ^ + (11 - 8c 3 )/8X. 

(1,2 ) 2 



(6.16) 



Here, ci and C3 are the c- values associated with the \k — h) and l/c = |) states. These 

2 2 1 Z' 1 Z' 

c-values are obtained as 



ci = |ci and cz = |ci. (6-17) 



Note that the P-baryon masses given by Eqs. (|6.17|) are quantitatively different from those 
of Table IV and furthermore that they violate the heavy quark symmetry; that is, there is 
no signal for any degenerate pairs of (1, \ ) and (1, | ). 

Then, what goes wrong in this straightforward extension? In the quantization procedure 
for obtaining Eq. (|1.2|) , only a single bound state is involved. Thus, to obtain the hyperfine 
energy of order 1/N C , it is enough to take the expectation value of H rot with respect to the 
unperturbed soliton- heavy meson bound state of order N®. However, in the heavy meson 
mass limit, all the heavy meson bound states come in degenerate doublets with grand spin 
k = ki± 1/2 (unless ke = 0) and consequently we have degenerate states up to order m^N®, 
for example, \i = l,j = |)) coming from k — | bound state and \i = l,j = |)) from k — | 
bound state. In evaluating the hyperfine energy of next order due to collective rotation, as 
is well-known in standard quantum mechanics, we have to diagonalize the energy matrices 
for the degenerate basis. As for our sample case, we obtain the energy matrices as 



11 c ( 8 9 FS\ 

£ (^ + ) = Mso1 + U ' B + 8J + 6Z (2^ -loj ' 

^>§ + ) = M °° l + J * + 8^ + f! (-2U ~-4 
which yield the P-baryon masses consistent with Eq.( |6.14| ) as 
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(6.18) 



m (1> i +) = M sol + u' B + (11 - 16ci)/8J, 

m 1+ = M sol + lo' b + (11 - 8ci)/8J, 

(1,a } (6.19) 

m (1) 3 +) = M soi + ^ + (11 - 8ci)/8Z, 

m (1) | +) = M soi + ^ + (11 + 8ci)/8J. 

It shows that the inclusion of the degenerate states into the quantization procedure is, 
thus, essential in restoring the heavy quark symmetry in the bound state approach. If we 
work with finite but sufficiently heavy meson masses, we have only approximate degenerate 
states of the grand spin k = fc^± |. Although the mass corrections due to the nondegenerate 
states are at most of order 1/N^, because of the small energy discrepancy of order 1/mg 
in the denominator, their couplings can compete with diagonal terms of order m^A" 1 and 
should be taken into account properly. 

VII. SUMMARY AND CONCLUSION 

In this paper, we have discussed the heavy quark symmetry in describing heavy baryons 
containing a single heavy quark or antiquark as bound states of the SU(2) soliton and 
heavy mesons. We have developed a consistent bound state approach so that the heavy 
quark symmetry is realized explicitly in the heavy baryon spectrum in the infinitely heavy 
mass limit. The resulting mass formula reads 



3 

m (o,D = M soi + UJb + ]pf> U* 

m (ij*) = M ™i + u B + — , cr 



82" 



2,2 (7.1) 

2 ' 2 / 



for the heavy baryons and 



3 
m(pj*) = M sol +u' B + —, (f 



rrnij*) 



8X 

M so1 + uj' b + ^, (^ = i ± ,f 



I ± 3±A 

O 1 O I 

(7.2) 



for the P-baryons. All the masses are consistent with the heavy quark symmetry and the 
genuine hyperfine splittings vanish. As for the heavy baryons, one can still apply the mass 
formula ( |1.2| ) to obtain their masses and the c-factor defined by Eq. (|1.3|) plays the role of 
the hyperfine constant. For the P-baryons, the degenerate states up to order mnN® requires 
diagonalizing the Hamiltonian with respect to the degenerate basis to obtain the hyperfine 
energy of next order in 1/N C . In this case, the nonvanishing c-factors defined by Eq.( |1.3|) 
should not be treated as the hyperfine constant. In restoring the heavy quark symmetry 
in the P-baryon spectrum, the couplings of the states made of the degenerate heavy meson 
bound states with the grand spins k = tt£ ± 1/2 are shown to play an important role. 

Except the ground state with i=0 and j 7r =| , the occurrence of the parity doublets in the 
spectrum (i.e. , two states with the same angular momentum but opposite parity occurring 
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at the same mass) is interesting. However, we are not in the position to conclude whether 
such a parity doubling has any physical importance as discussed in the context of chiral 
symmetry [32 1 and Regge-pole theory [33] or just an artifact from our approximation on 



the radial functions. To extract more decisive conclusion, we should work with finite heavy 
mesons incorporating the kinetic terms MM and see whether the parity doubling occurs in 



the heavy quark limit. Iachello [45] has reported that a similar parity doubling in the excited 
baryon spectra is occurred in the baglike models and stringlike models and it was analyzed 
as a consequence of the geometric structure of baryons. 

We have corrected a mistake committed in applying the traditional bound state approach 
of evaluating the expectation value of the operator 2 to the heavy baryons. It has been 
approximated simply by the square of the expectation value of as 

((iJe\e%j e ))Kc%(k e + l). (7.3) 

In our approach, it can be exactly obtained as 

(^,M® 2 \hJi)) = i (7.4) 

In case of heavy baryons with c = 0, the correction is an overall shift of the heavy baryon 
masses by an amount 3/8X ~ 60 MeV. 

Furthermore, the simple structure of the model enables us to illustrate explicitly how the 
spin-grandspin transmutation occurs in the bound state approach. In the isospin co-moving 
frame, the total spin of the soliton-heavy meson bound system can be obtained as 

J = R + K M , (7.5) 

— * — * 

with the rotor spin R and the grand spin K^f of the heavy meson fields in the isospin co- 
moving frame. The latter plays the role of the heavy-meson spin; that is, the isospin of the 
heavy mesons is transmuted to the part of their spin. 

We have worked with infinitely heavy mesons (and soliton). It provides a useful instruc- 
tion to the bound state approaches with finite-mass heavy mesons so that it has a correct 
heavy quark limit; that is, one has to include the nearly doubly degenerate heavy meson 
bound states of grand spin k = kg ± 1/2 into the quantization procedure. 
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APPENDIX A: SPIN AND ISOSPIN OPERATORS AND THEIR EIGENSTATES 

The invariance of the Lagrangian density ( |2.2|) under the infinitesimal Lorentz transfor- 
mation 
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x'» = x" + e(X, with e"" = -c"" 



-i/ - ' > 



$ _► $'(a/) = $( z ), (Al) 

d>*(x) - $' a *(x') = ^(S^J^x), with (S^ = 3%-// a 

defines conserved angular momentum operators 
J^le^f^rMojk, 

M ojk = f x j<pok _ x k V o^ + (n* m (s fcj ') m „$* nt + ^* m (s kj ) mn u* n ^ 



(A2) 



where P Miy is the canonical energy-momentum tensor and IT n (= 9£f ree /9$*) is the momen- 
tum conjugate to the field <3>* n . Here, the indices run from 1 to 3. The first part corresponds 
to the orbital angular momentum, L, and the second part correspond to the spin angular 
momentum, S. In the heavy meson mass limit (m$ $* — > oo), substitution of Eq. (|2.3j a) with 
v = (1,0) leads to the spin operator of the heavy meson fields as 

S = % [d 3 r $* v x $$ (A3) 



as the leading order term in the meson masses. From now on, we will work in the rest frame 
of the heavy mesons. In terms of the 4x4 matrix field H(x), it can be simply rewritten as 

S=- fd 3 rTr([y,H}H), (A4) 



which implies that the corresponding quantum mechanical spin operators in the 4x4 matrix 
representation is the Dirac spin matrices acting on the wavefunction H n (x) as 

S{H n }=[\a,H n }. (A5) 

The minus sign in Eq. (|A4|) is due to the normalization convention of the field H(x) and 
H{x). Note that the meson number operator is given by 

N = - fd 3 rTi(HH). (A6) 

One can easily find the eigenstates of the spin operators as 

( S =0; S3 =0) = [^(l + 7°)]75, 

(^=1^3=0) = [^(l + 7 )]7 3 , (A7) 

(s=l;s 3 =±l) = [^(1 + 7°)H5§(7 1 ± *7 2 )], 

which are normalized as 



Tr{(s;s s )(s>;s> 3 )} = -8 S , S 5 S ,,. (A8) 



Furthermore, the invariance of the Lagrangian (|2.4j) under the heavy-quark spin rotation 
3|b) leads to conserved heavy-quark spin operators Sq as 

24 



S Q = - d 3 r Ti{l& HH), (A9) 



that is, the action of the heavy-quark spin operators on the wavefunction is the multiplication 
of the Dirac spin matrices to its left-hand-side: 

S Q {H n } = \oE n . (A10) 

Since the heavy-quark spin decouples in the heavy meson mass limit, it is convenient to 
introduce the "light-quark" spin operators for the light degrees of freedom in the heavy 
mesons: 

S e (= S-S Q ) = + jd 3 rTr{H\aH), 
Si{H n } = H n {-\a). 

The eigenstates of Sq and Se can be explicitly obtained as 

^(+111 + |)q = (s=l;s 3 =+l), 

<(+ill - §)q = 7l{(s=l;s3=0) + ( S =0; S3 =0)}, 

.(-III + \)q = j= 2 {(s=l;s 3 =0) - (s=0;s 3 =0)}, 

On the other hand, the isospin operators associated with the invariance of the Lagrangian 
Q2.15 ) under the isovector transformations ( |2.8|) and (|2.12|) (with L=R=U) are 



I = I M + h, (A13) 

with I M the isospin operator of the Goldstone boson fields 

f 2 



I M = i Jd 3 r^ Tr{ir(St,9 S + £<%£+)} + • ■ ■ , (|Al3| a) 

and Ih those of the heavy meson fields interacting with Goldstone pions 

h = +iJd 3 rTr{men + ^)H}. O) 

Especially, the isospin operator of the free heavy meson fields is 

d 3 rTi{H(-±T)H}. (A14) 



Thus, the quantum mechanical isospin operator for the wavefunction is the 2x2 Pauli 
matrices (— |t) acting on the right-hand-side of the anti-isodoublet H n (x): 

I h {H n {x)} = H n {x){-\f). (A15) 

Explicitly, their eigenstates can be written as 

0+1 = (0,-1), and 0_i=(l,O). (A16) 
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APPENDIX B: KVBASIS 



Here, we present the explicit expressions of the the angular part of the wavef unctions, 
which are the eigenstates of Kg, Kg jZ and Sq, Sq jZ . They can be obtained by linear combina- 
tions of direct products of four angular momentum bases; Y^ m (f)(the eigenstate of L 2 , L z ), 

i (the eigenstate of I\, Ih, z ), ^(±§| (the eigenstate of Sf, Sg >z ), and | ± |)g(the eigenstate 






of Sq, Sq :Z ). For a given quantum numbers, ke,k 3 , we have four states K k j k , i=l,2,3,4 
according to the numbers of the different combination of Ih and Sg. 

We first combine the spherical harmonics Yp m (r), and the isospin basis ± i to obtain 

y\x 3 (r), the eigenstates of A 2 and A 3 (A = L + I h ): 



,-! / A-A 3 + l ^ 7 / A + A 3 + l . 

A ' Vi "V 2(A+1) W±W+± + V 2 (A+1) ' ,AH -- W '"-+ " " 



^ = \/-^^A,-i(r)^ + i + \/-^^A, + i(r)^-i, A = £ + i 

(Bl) 



2 2 



It provides a convenient basis in evaluating the expectation values of the operator including 
(r-r), since y^ satisfy a useful identity 

y^-?) = y$ 3 - (B2) 

Next, the eigenstates /Q^ fc3SQ is obtained by combining ^W and £(±||| ± §)q : 
(i)i=l; A=£+|, A;,=A-i=£ 



(1) _ kg-k 3 + l (+) x kg + k 3 + l (+) x 

***»•* " "V 2(fo + l) JW-^+sIWq + V 2(^ + 1) ^-^"IMc 

(ii) i=2; A=£ - i A; £ =A + i=£ 



r( 2 ) / ^ + fc 3 ^,(-) /,i M x kg-h ( _ } x 

(iii) j=3; A=£ - |, fc £ =A - |=£ - 1 



(B3) 



(3) _ kg-k 3 + l (_) x kg + h + 1 (_) x 

^fc/fcs-o- V 2(fc, + l) ^3-^1+211^0 + y 2(fc, + l) ^3+i^ 2ll s QJ 
(iv) i=4; A=£ + i, fc,=A + |=^ + 1 



r (4) _ 1 1 / ^ + fc 3 v (+) /,i|i_x 1 i / ^- fc 3 v (+) / in 



They are normalized as 

One can easily check that K/ k ^ kzs and IC kekas have the parity tc = —(—l) ke and the other 
two have the parity n = (— l) ke and they are related to each other as 
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K. 



(0 



ir-f) = K 



(i+2) 



'M*«qV ■' ; - '^ktk s s Q - (i-h 2 ) (B5) 

These .fQ-bases evaluate the matrix elements .My defined by Eq. (]3.10D as 

L 



M 



M 



jggCg) 

2ife, + l 



yF'(O) 



2 



1 o 1 

1, 



(fori,j=l,2) 



2k e + 3 -4Jk e (k e + 1) 



/ 



\-^h{ki + 1) 2fo - 1 



(B6) 



(for i,j=3,4) 



independent of k% and the heavy quark spin Sq. 

In Sec. V, we work with the expectation values of the operator B(oo) with respect to 

— * 

the single-particle Fock states \n) for which we need to evaluate the expectation values of 9 



with respect to the .fQ-basis. Wigner-Eckart theorem enables us to express them as EB 



,<*&|S«(°°)IM*>< = - /dniKK® r -rilrYr-fm. } = foJ^pS> (Afl|3||fr)< 



2fc; + i 



(B7) 



with the "reduced matrix elements" : 



h)i 



x {h 


| <S> || 


2 (ke 


|9|| 


s(kt 


|e|| 


s( k t 


|9|| 


i{kt 


l©ll 



1 k e (2kt + l) 



h + 1 



1 (k e + l)(2k e + 1) 

2V fc, 

1 (2^ + 3)y/k~i 

2 {h + i){2k e + iy 

1 



4 ~ V2h + T\ 

1 (2kt - l)^/hTl 



k t {2k t + 1) 



(B8) 



^-llieHA;^ 
x (h- 1||9||^) 4 

2 (h + 1||9||^) 3 



'(A* + l)(2A;i-l) 





2A; £ + 1 


1 


(2ki - 1) 


2 


\ 


k e (2kt + lY 


1 


(2k e + 3) 


A 


( 


k e + l){2k e +l) 



j(^ + l||e||A;. 



£) A 



t h(2k e + 3) 
2k e + l 



and others are zero. As far as the single-particle Fock states of the same kg are concerned, 
we can rewrite Eg. (|B7|) in a more convenient form as 
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a {ki,m'\Q{oo)\ke,m) b = -c ab (k e ,m'\K e \k E ,m), (B9) 

where (kg, mf\Kg\ke, m) denotes the expectation value of the operator Kg with respect to its 
eigenstates. The multiplication coefficients analogous to the Lande's ^-factor can be read 
off from Eq. (|B~7|) as 



C 3S — I 



2k e + 3 2h - 1 



33 ~~ t o/l , ,\/ol i i\> u 44 



2(k e + l)(2k e + l) 44 2kt(2kt+l) 
c, 
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'11 — o/; i i \' u 22 " 



fc £ (fc £ + f ) 2k i + 1 
f f 

2(A*Tl) ' ° 22 ~ + 2fc? 



and others zero. With respect to the states \ki)±, we have 

f 



0, c ++ = 4 



2 JA*(fo + f) 
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2kt(k e + iy /IP[TWX 

1 v ; (|B|b) 
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FIGURES 



FIG. 1. Bound states of (a) heavy mesons and (b) antiflavored heavy mesons. 
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FIG. 2. Heavy baryon spectrum. 
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FIG. 3. P-baryon spectrum. 
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TABLES 



TABLE I. Heavy-meson eigenstates with £ = and 1. 



{k e ,k 3 ,TT,s Q } 


£ 


eigenfunct. 


4ff 


fc 


n) 


{0,0, -,s Q } 


+bF'(0) 


f(r)JC^s Q 


1 


l 

2 


|0)i 


{l,m,-,s Q } 


-lgF'(0) 


f(r)1ctl Q 


1 


1 3 

2' 2 


|1>- 


{l,m,-,s Q } 


+ l 2 9F'(0) 


/( r )^lros (3 


1 


1 3 

2' 2 


|1)+ 


{o,o,+, Sq } 


-lgF'(0) 


f(r)K§L Q 





1 

2 


|0>3 


{l,m,+,s Q } 


HgF'iO) 


/(o4L Q 
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1 3 

2' 2 


|l)l 


{l,m,+,s Q } 


+bF'(0) 


f(r)K?L Q 





1 3 

2' 2 


|1>2 


{2,m,+,s Q } 


-lgF'(0) 


f(r)JC { iJ SQ 


2 


3 5 
2' 2 


|2>- 


{2,m,+,s Q } 


+yF'(o) 


f(r)4~£s Q 


2 


3 5 

2' 2 


|2>+ 



TABLE II. 



states for heavy baryons. 



i 


tf 


n) 


l*,tf» a 


£ 


3 







0+ 


|0> 


|o,o+» 


-lgF'(0) 


i 

2 


Aq 


1 


i+ 


|0) 3 
|2>- 


lLi + )) 2 


-lgF'(0) 
-lgF'(0) 


1 3 

2' 2 


S <2' S Q 





l- 


|1>- 


|o,i-» 


-| 5 F'(0) 


1 3 

2' 2 




1 


o- 


|1>- 


IL0-)) 


-l^'(o) 


1 
2 




1 


1- 


|1>- 


|l,l-» 


-| ff F'(0) 


1 3 

2' 2 
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TABLE III. \i,jj)) states for the P-baryons. 



i 


ft 


\n) 


ri Jf //i 


£ 


3 





0" 


|0)i 


|o,o-» 


-bF'(0) 


i 

2 





1- 


|1>+ 


|o, i-» 


-\gF'{0) 


1 3 

2' 2 


1 


o- 


|1)+ 


|o, i-» 


-lgF'(0) 


1 
2 


1 


1- 


|0)i 

11)+ 

|2)i 

|2>2 


11, 1-»1 

|1,1 _ »2 
|1,1 _ »3 

|1,1~»4 


-bF'io) 


1 3 

2' 2 





1+ 


|l)i 

|1>2 


|0,l + ))i 

|o,i + )) 2 


-lgF'(0) 


1 3 

2' 2 


1 


0+ 


|1>2 


|l,0+))i 

|l,0+» 2 


-bF'(0) 


1 

2 


1 


1+ 


|1>2 

|2>+ 


IM + »1 

|M + »2 
|l,l + »3 


-\gF\Q) 


1 3 

2' 2 



TABLE IV. Positive and Negative Parity P-baryon (P) masses (in MeV). 



i 


3t 


j* 


Mass Formula 


m P 

c 


m p 


b.e.* 





cr 


i- 

2 


M sol +J B + 3/8J 


2704 


6042 


210 





i- 


1- 3- 

2 ' 2 


M soi + <J B + 3/8J 


2704 


6042 


210 


1 


i - 


1- 3~ 

2 ' 2 


M soi + u4 + 7/8J 


2802 


6140 


112 





i+ 


1+ 3 + 
2 ' 2 


M soi + J B + 3/8J 


2704 


6042 


210 


1 


0+ 


1 + 
2 


M soi + u^ + 7/82 


2802 


6140 


112 


1 


1+ 


1+ 3 + 
2 ' 2 


M soZ + lo' b + 7/8J 


2802 


6140 


112 



* Binding energy below M^. 
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